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As the result of recent interest in the possible weight saving obtained by using high contact ratio 
gear (HCRG) meshes, a computer program was developed using a time history, interactive, closed 
form solution for the dynamic tooth loads for both low and high contact ratio spur gears (ref. 1). 
Because the magnitude and variation of the tooth pair compliance with load position can affect the 
dynamics and loading significantly, and because the tooth root stressing per load varies significantly 
with load position, it was desirable to include in the dynamic gear program preprocessors and 
postprocessors, respectively, for calculating these two important factors. Since the tooth forms for 
HCRG deviate appreciably from the tooth forms for conventional low contact ratio gears (LCRG), 
improved and simplified methods were developed for calculating the compliance and stress sensitivity 
for three involute tooth forms as a function of load position — a standard LCRG tooth with no 
undercut and two HCRG teeth with different forms of undercut (see fig. 1). 

The method developed for calculating the compliance of spur gear teeth follows to a great extent 
that developed by Weber (ref. 2) and includes three factors: (1) the basic deflection of the tooth as a 
beam, (2) the deflection of the tooth caused by the fillet and foundation flexibility (ref. 3), and (3) the 
local deflection caused by the contact between the two teeth. The principal improvement in the 
compliance analysis was for the fillet/foundation deflection, which was found to be defined by 
different fillet angles for various load positions; whereas previous studies assume a given fillet angle 
of about 75°. The resulting compliance analysis was evaluated by applying it to available test, finite- 
element, and analytic -transformation results and was found to give compliance results which agreed 
well with measurement and “exact” analyses. 

The method developed for calculating the stress sensitivity is an improved and simplified version 
of the Heywood analysis (ref. 4). Because the tooth forms for HCRG deviate appreciably from 
conventional LCRG tooth forms, the stress sensitivity analysis had to include most of the factors 
affecting the stressing. The improved analysis does not assume the peak stress occurs at 30° from the 
base of the fillet, as the Heywood analysis does, but allows it to be a function of load position. The 
sensitivity analysis was found to give results which agreed well with the Heywood (ref. 4) and Kelly 
and Pederson methods (ref. 13), which, in turn, have been found to correlate well with test results. 
Evaluation of the modified Heywood stress sensitivity analysis showed that its results compared very 
well with available test, finite element, and analytic transformation results. The differences between 
the improved Heywood type stress sensitivity analysis and the conventional Lewis and AGMA 
analyses were found to be significant for low and high contact ratio gears, being about 25 percent 
higher for the LCRG and about 15 percent higher for the HCRG. 


Gear Tooth Compliance 

Derivation of Compliance Formula. 

The determination of the compliance of gear teeth is considerably more difficult than 
determining the stress sensitivity because it is an integral function of the entire loaded tooth, rather 
than just a function of the section properties at the peak stress point and the load location. In 
addition, because of the stubbiness of the teeth, the foundation and shear effects are important. The 
method developed herein parallels to a great extent Weber’s work (ref. 2) but using O’Donnell’s 
foundation factors (refs. 3 and 5). The total compliance or flexibility of a gear tooth at the point of 
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Figure 1. - Gear tooth geometries. 



load, yj, is made up of three deflections: (1) the basic tooth as a cantilever beam, y^-, (2) the fillet and 
foundation, yp\ and (3) the local contact and compression, yp. The compliance, C, of a gear tooth 
pair, 1 and 2, is their combined deflection per unit of load at the contact position or 

c= (>n +y-nVL = [{yg^ +yBi) 

+ (yn+yn)+(yL)VL (i) 

where L is the applied load and the deflections y are in the direction of the load. The derivation of the 
three factors in equation (1) follows. 

In order to determine the beam compliance of a gear tooth its geometry must be defined over its 
entire load length. Three different involute tooth forms are considered in this study, which depend on 
the contact ratio, number of teeth and pitch, and manufacturing process. The three involute tooth 
forms are depicted in figure 1 , form S being for standard LCRG teeth with no necking and forms A 
and B being for HCRG and small radium LCRG teeth. For tooth form S the tangency of the fillet 
radius with the involute form occurs outboard of the base circle. For forms A and B the tangency 
with the tooth form occurs inboard of the base circle. The tooth form inboard of the base circle is 
machined along a line tangent to the involute surface at the base circle for form A, whereas it is 
plunged ground along a line parallel to the fillet radius centerline for form B. The appendix gives the 
expressions for defining the tooth geometry for three forms of teeth. The applicable tooth geometry, 
form S, A, or B in figure 1, is determined by the expressions 


FormS:(/?3 + 2r/?,) a 

Form A or B: (R} + Ir R,) < Rl (2) 


Whether form A or B applies is defined by the manufacturing process used, that is, fillet region made 
tangent or plunge ground, respectively. Having analytically defined the geometry of the gear tooth, 
the compliance of the tooth can be determined. 

The deflection and, therefore, compliance of a gear tooth over its beam portion is easily 
obtained using elementary strength of materials. Referring to figure 2(a), the total bending and shear 
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deflection in the direction of and at the applied tooth load, L, which is at radius or position S 
along the line of action, can be expressed in integral form as 


L Cos^ <l>i f^o 
“ £ Jo 




(3) 


where / is the section modulus of the tooth as a function of^=xor« = z;« = (/-0 = z = (/-jc)andis 
the distance along the centerline of the tooth from the load position; a value of 1.2 has been assumed 
for the shear factor based on a rectangular tooth; G = £'/2(l + ^); and A is the cross-sectional area as 
a function of x or z. The deflection of the basic tooth can also be defined in a summation expression 
rather than integral form (see fig. 2(b)), which would be more applicable for a hand calculation or a 
programmable calculator. In this case the tooth beam deflection at and in the direction of the load is 


Xb = 


L Cos^<t>L y 




(//-/, 5 ,+ ^ 6 ,^) 


(2.4(1 +M) + Tan^0i) 

A, 


(4) 


wher^l/7j- = (!//,• -I- 1//,+ 1)/2 and (1/A,+ 1/A,+ |)/2. Using these inverse forms for the values of /,■ 
and Ai improves the accuracy for a small number of elements. In equation (4) Ij = (i-Xi) and 
dj = (Xi + 1 -x,)(see the appendix). Both approaches for beam flexibility assume a narrow tooth width, 
W. For wide teeth where W/hp>5, the flexibility is decreased by the anticlastic effect, so that the 
values of I in equations (3) and (4) should be divided by (1 -yfl)- 

Because of the fillet and the flexibility of the material to which the tooth is attached, additional 
deflection will occur at the load (refs. 2, 3, and 5). This fillet and foundation deflection in the 
direction of the load, yp, is a function of the fillet geometry and the load position and direction and is 
determined by the effective fillet length or angle 7 /rfor which the maximum deflection or work occurs 
at the load. In other references this effective fillet angle has been assumed constant at about 70 (ref. 
3) or neglected (ref. 2). 

Based on figure 3, O’Donnell (refs. 3 and 5) shows that the deflection at and in the direction of 
the tooth load due to the foundation effects, ypp, for plane stress is 


L Cos^ r 16.67 ( If y 

WE L r \hf) 


+ 2 ( 1 - m )(^) + 1 . 534(1 + 


2.4(1 +/i)/J 


(5) 



Figure 3. - Fillet and foundation compliance of a 
gear tooth. 
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For wide teeth the expression for plane strain is used or 


yfF = 


L Cos^ ()>i 
~WE 






( 6 ) 


The O’Donnell coefficients in equations (5) and (6) differ slightly from those given by Weber (ref. 2). 
The first term in the brackets is the deflection at L due to the rotation caused by the moment at hp. 
The second term is the sum of the deflections at L due to the displacement dX hp caused by the 
moment at hpand the rotation at /i/r caused by the shear force at hp. The first part of the third term is 
the displacement at L due to the shear force at hp, based on the assumption that the effective depth 
for determining this deflection is 2Vi times the tooth thickness (ref. 2). The second part of the third 
term is the deflection at L due to the normal component of the load assuming the same relationship 
holds as indicated by the beam deflection equation (see eqs. (3) and (4)). The question now follows: 
What values of Ip and. /i/r should be used in equations (5) and (6), that is, what is the effective fillet 
length or angle -ypl 

Referring to figure 3, the deflection at and in the direction of the load due to the flexibility of the 
fillet and foundation, yp, the shaded region, is obtained from the summation of the fillet beam 
deflection, ypp, using equations (3) and (4), and the foundation deflection ypp, using equation (5) or 
(6), that is, 

yF=yFF+yFB O) 

where //r=r+r(Sin 7/7-Sin 7), hp= h + 2r{Cos y- Cos yp), /, = r+r(Sin 7,-Sin 7), and hi=h 
+ 2r(Cos 7 -Cos 7,). The value of 7/7 is the one that maximizes the value of ypor yp, which can be 
easily done as one integrates or sums up the deflection of the tooth starting at the beginning of the 
fillet or at the load position. 

The significance of maximizing the deflection due to fillet effects can best be realized by 
expressing in nondimensional form the deflection at the load due to the fillet and foundation 
flexibility^ If this is done, it is found that the effective fillet angje, yp, is a function of four quantities 
-T/Ti, r/h, 7, and Figure 4 shows the effect of r/h and l/h on 7/7 for typical values of 4>i = 20° 
and 7 = 0. Apparent from these curves is that assuming a value of 75° for 7/7 (ref. 3) applies only for 
small values of J/h and r/h and that a more universal value would be 55° to 60°. The effects of 7 and 
4>l on yp were studied and found to be small. Figure 5 presents the corresponding fillet and 
foundation flexibility results. Apparent^ from these curves is that assuming a value of 75° for yp 
results in progressively greater error as l/h and r/h increase (HCRG). However, for typical values of 
J/h and r/Ji, particularly for LCRG, the likely difference will be less than 10 percent from that given 
by solving for yp. 

The local compliance, yp, consists of the Hertz or line contact deflection plus the compression of 
each tooth between the point of contact and the tooth centerline. Figure 6 gives the nomenclature for 
the parameters that determine the local deformation. Three viable approaches were reviewed with 
regard to calculating the local compliance: (1) an approximate Hertzian and compression approach 
originally used at Hamilton Standard, (2) a semi-empirical approach developed by Palmgren (ref. 6), 
and (3) a closed form approach developed by Weber (ref. 2). 

The approximate approach assumes that the true Hertz contact deformation, >>//, must be 
doubled to account for the cupping action discussed in reference 7, that is. 




2 



4/ ( \-n} 

tW V £, £2 / 


( 8 ) 


The local compression of each tooth between the contact point and tooth centerline is approximated 
by assummg the load spreads at a 45° angle, that is from 2b on the top to 2(b -I- h) at the bottom over 
a length h. Assuming the compression deflection is given by the average width, we find 


" trlV £, /\2b + h,/ V £, /\2b+hJi 


(9) 
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Figure 4. - Fillet angle for maximum flexibility 
versus load point from fillet. 



Figure 5. - Flexibility of fillet and foundation 
versus load position from fillet. 
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Figure 6. - Nomenclature for local compliance. 






where 





(10) 


Assuming hi/b and h 2 /b are both greater than 1, the local compression of the two teeth reduces to 


yc 



)] 


so that the overall local compression of the two teeth is 


( 11 ) 


vl =yH +yc 





which for £’1 =£2 and ii\- 1^2 reduces to the simple expression 


( 12 ) 




»3.57 


4(1 -M^) 
tE 


L 

W 


«4.55(1-m^) 


L 

WE 


(13) 


Although equations (12) and (13) show_y/^ increases linearly with L/W, in actuality it increases at a 
slightly lower rate because the effect of the contact width was neglected. This nonlinear effect can be 
included if desired. 

Reference 6 presents an expression for the local deformation based on the semi-empircal 
equation developed by Palmgren for contacting cylinders in roller bearings. Because only half a tooth 
thickness is being compressed instead of the entire roller, the local tooth deflection was assumed to be 
half that for the roller bearing or 


1.275 L ’ 


4(1-m^)£ 

tE.^W 


( ) 


( 14 ) 


where l/£i 2 = (l/£'i + l/£’2)/2. Weber in reference 2 developed an expression specifically for the 
local deformation of two gear teeth. In order to obtain a closed form solution, he assumed small 
deformations so that just the first two terms of the binomial expansion of the deformation needed to 
be used, that is. 



where b is given by equation (10). If E\ =£2 = equation (15) reduces to 


yi = 


4(1- 




(16) 


All of the expressions for the local deformation are nonlinear with load because of the Hertz half 
contact width b. For simplicity the nonlinearity can be circumvented by assuming a value of b based 
on the maximum nominal value of tooth load during the mesh, that is, T/CRmin-Cos <!>pRp, or 
assuming a maximum value of Hertz stress. This latter option should be used after the initial dynamic 
gear tooth analysis is run and a first approximation of the maximum Hertz stress is obtained 
including dynamic effects. 

In order to decide which of the three local compliance methods should be adopted, they were 
evaluated for several cases. The first case was for a pair of steel, standard 2-in. pitch radius, 8 pitch, 
0.375-in. wide, 20° pressure angle gears in contact at the pitch radius (/i = 0.0998 in.) under load 
levels of L/W= 1279 and 2875 Ib/in {q= 140 000 and 210 000 psi), respectively. The second case is 
from reference 2 and was for the meshing of an 18 tooth, steel gear with a rack under three loading 
conditions, where =0.80 m, ^2 = 0.84 m, /-j = 3.1 m, and /- 2 = <». Here, m is the tooth module-pitch 
diameter in millimeters divided by the number of teeth. The local deflections predicted by the three 
methods for these two cases are summarized in table 1. For the first case, the empirical method gives 
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Table 1 Comparison of Methods for Calculating Local Tooth Compliance: Approximate; Empirical (Palmgren); 

Closed Form (Weber) 


Factor 

Equation 

Gear Caaa 1 

Gear Caaa 2 

16 Taath/16 Taath 
8 PiKh. W • 0.375 in.. Rp - 2.0" 

E - 30 X 10* pai; <pp • 20°; h - 0.0998" 

18 Taath/ o^Tmth 

• 0.80 m.^2 * 0.84m; r^ ■ 3.1m; f 2 ■ 

E • 2.15 X 10* Kg/cm^, m ■ Module in mm 

Load/Witfth 

L/W 

1279 Ib/in. 

2875 Ib/in. 

16.16 Kgnt/cm^ 

23.00 Kgm/cm^ 

32.33 Kgm/cm^ 

Hartx Siraaa 

q 

140 ksi 

210 ksi 

1400 Kg/cn<2 

2000 Kg/cm^ 

2600 Kg/em^ 

Half Contact Width 

b 

0.00581 in. 

0.00871 in. 

0.00735 m 

0.01050 m 

0.0147 m 

Local 

Appron-Linaar 

(131 

3.57 

3.57 

3.57 

3.67 

3.57 

Complianca 

Approx-Non Lirtaar 

(9) 

3.41 

3.34 

3.52 

3.48 

3.44 


Empirical 

(14) 

2.50 

2.31 

3.86* 

3.73* 

3.59* 

V ffEW 

Clotad Form 

(16) 

3.32 

2.91 

5.19 

4.83 

4.46 

4L (1 









AMumtd Equal to 5m 


results about 80 percent of those by Weber’s method, probably because it includes a width effect, 
whereas the approximate method gives 10 to 20 percent high values, particularly for the higher Hertz 
stressing. For the second case there is considerable difference in the local deflections predicted by the 
three methods for these relatively low Hertz stress cases. Weber’s expression gives values about 20 to 
25 percent greater than the other two. Because the empirical method is width dependent and appears 
to give values somewhat low and because the approximate method does not appear to include the 
nonlinear effects properly, it was decided to use the closed form expression of Weber’s equation (15), 
for the local compliance. 

The overall compliance, C, of the tooth pair is obtained by adding the local gear tooth 
compliance defined by equation (15) with the gross tooth compliances for the two gear teeth 
(>'a+>'/r)/Z, given by equations (3) or (4) and (7). (See eq. (1).) In order to obtain the necessary 
compliance information for the dynamic gear analysis of reference 1 , the tooth pair compliance must 
be determined for at least five different load contact positions, S, along the line of action, so that the 
compliance can be expressed in the required five term power series 

C=C„[l+>l(S/SJ 

+ fl(5/5j2 + c(S/Sj3+Z>(5/SJ“] (17) 

Here, Sg is the reference distance along the line of action corresponding to the tip radius of the gear, 
Cg is the tooth pair compliance for a load at the pitch radius, and A, B, C, and D are curve fitting 
coefficients. 


Evaluation of Compliance Formula 

The cantilever beam test and finite elements results of reference 8 provided a means for checking 
the basic fillet and foundation formulation, because of the three deflection measurements along the 
beam, the simplicity of the beam itself, and the three different fillet sizes and beam thicknesses. (See 
fig. 7.) A summary of the comparison of the finite element, O’Donnell, and formulas (4) and (7) 
theoretical results with test is given in table 2. Although the latter formulas are based on O’Donnell’s 
foundation analysis, the effective fillet length is solved for based on maximum energy, resulting in a 
fillet length equivalent to an angle 7 fOf about 32° to 34°, depending on the case. Table 2 shows that 
all three theoretical methods do quite well at predicting the load deflection; however, the O’Donnell 
approach appears to be inaccurate for small values of r/ hg and Ig/ hg and for load positions close to 
the support juncture or small x. The error for the formulas herein tends to be more uniform for the 
various tooth geometries as does the finite element results. Correcting the test results for the Poisson 
swelling (see ref. 8) makes the error variation with position jc nearly uniform at about +1 ±3 percent 
for the results based on formulas (4) and (7). However, because of the large load moment arm, the 
concept of variable effective fillet length developed herein could not be fully confirmed. 
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JiH 

Fillet Radius — 

R - In. 

Beam Thickness - 

- H - In. 

R/H = 

R/H = 

R/H = 

R/H = 

R/H = 

R/H - 


0.12 

0.31 

0.5 

0.12 

0.3 

0.5 

S 

0.188 

0.500 

0.812 

1.60 

1.60 

1.60 

8 

0.125 

0.313 

0.500 

1.00 

1.00 

1.00 

10 

0.094 

0.250 

0.375 

0.80 

0.80 

0.80 


Figure 7. - Geometry of cantilever test specimen 
(ref. 8) . 


Table 2 Percent Error in Deflection Compared to Experimental Data; Cantilever Beam: Load at / = 8", 
Deflections at X, Thickness H, Radius R; reference [8] and Fig. 9 


Configuration 

X - 1.5" 

X - 3 .0 " 

X - 5.0" 






Equations 

1 


Equations 



Equations 

R/H 

t/H 

Equation 7 

O'Don nail 

F.E. 

4&7 

O'Donnall 

F.E. 

4&7 

O’Donnall 

F.E. 

4&7 


5 

34 

+6 

S 

+2 

+7 

-4 

+1 

+1 

-6 

+4 

0.50 

8 

34 

+5 

-3 

0 

0 

-4 

0 

+4 

+1 

0 


10 

34 

+3 

-3 

-2 

+7 

+4 

0 

+4 

+2 

-4 


5 

32 

-5 

-4 

1 

-3 

-3 

-1 

0 

0 

3 

0.31 

8 

32 

-3 

+1 

-4 

+2 

+4 

+4 

-2 

-1 

-2 


10 

32 

-4 

0 

-3 

-2 

0 

_4 

+1 

+2 

+1 


5 

32 

-11 

-6 

+3 

-6 

-3 

+5 

0 

+1 

+3 

0.12 

8 

32 

-8 

-3 

0 

-5 

-2 

+3 

0 

+2 

+1 


10 

32 

-7 

-2 

-2 

-4 

-2 

-1 

+1 

+3 

0 

Avaraga 

it Variation 

-3 ±8 

-3 ±5 

-0.6 ±3 

-0.5 ±6 

-1 ±5 

0.8 ±4 

+1 ±3 

0.5 ±3 

1 ±3 


In order to evaluate the compliance method herein for predicting the effects of axial 
compression, partial fillet, and beam stiffness present in actual spur gear teeth, it was applied to a 
number of the tooth configurations analyzed by complex, analytic transformation (refs. 9 and 10) 
and by finite element analysis (ref. 11). Reference 9 analyzed the beam deflection of the tooth 
centerline at various load locations for standard 20° pressure angle, full depth gears with 20 to 80 
teeth and compared the results to those predicted by Weber (ref. 2). The results, which are for plane 
strain, are reproduced in figure 8 along with the plane strain results predicted by the beam 
compliance method herein. The shape of the curves with load position along the line of contact, S/P, 
and the relative effect of number of teeth are similar for all three methods; however, there is quite a 
difference in their normalized compliance, yElV/L. In general, the transformation compliances (ref. 
9) are about 50 percent of those by Weber and about 60 percent of those predicted by formulas (4) 
and (7). 

Very important to the accuracy of the transformation method is the reference point used to 
define the deflection. If this is chosen close to the base of the tooth, low relative deflections will 
result. Reference 9 does not define the deflection location. The importance of the reference location 
is illustrated by the finite element results of reference 1 1 for 17 and 45 tooth gears, which are plotted 
in figure 8. A review of the investigatory results in reference 1 1 shows that by restricting the length of 
the foundation of the finite element model to that between the fillet centerlines, the tooth deflection 
at the tip was decreased at least 10 percent. If the results of reference 1 1 are corrected accordingly, the 
finite element compliance results for the 45 and 17 tooth gears would be about 85 and 95 percent. 
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Flexibility - YEW/L 


respectively, of those predicted by the method herein. Also, the shapes of the compliance curves are 
more consistent with the other results, for the added correlation increases as S/P becomes more 
negative. Figure 9 gives the derived values of effective fillet angle, 7 /r, for the various load positions 
and gear teeth which will result in the maximum deflection or work done by the applied load. For 
load positions inboard of the pitch circle, the assumed value of 75° by O’Donnell is a good 
approximation; however, for load positions towards the tooth tip, particularly for HCRG, the 
effective angle should be considerably less resulting in greater flexibility. 

As a further check on the accuracy of the beam and fillet compliance formulas (4) and (7), they 
were used to predict some of the tooth pair beam stiffnesses presented in reference 10 , which used the 
same analytic transformation method as reference 9. Because the results included the local effects 
based on Weber (ref. 2), they were corrected for this effect. Table 3 summarizes the results and the 
corresponding results from reference 9 and formulas (4) and (7). Apparent is the very good 
correlation of the results of this work and reference 10, whereas those of reference 9 (fig. 8 ) appear 
low. The differences between the results of references 9 and 10 are strange because they are based on 
the same transformation equations; however, the reference points for defining the tooth deformation 
might be different. In general, the results using the compliance beam formulas herein and those of 
reference 10 agree within 20 percent, and for the tooth pair they agree within 1 percent. 

Reference 12 presents the measured stiffness of a pair of standard form, steel gear teeth with the 
load at the pitch radius. The 20° pressure angle gears were 14 -in. wide and had a diametral pitch of 3, 
one with 27 teeth and the other with 18 teeth. The measured stiffnesses at 1275 and 1492 lb tooth 
loads were 1.28 x 10® and 1.49 X 10^ Ib/in, respectively. The theoretical method herein gives tooth 
stiffnesses of 1 .37 and 1.38 x 106 ib/in, respectively. The local flexibility represented about 25 percent 
of the overall flexibility. These results compare very favorably with the measured values, being only 
about 1 percent greater than their average value. This, coupled with the good correlation with the 
tooth pair results of reference 10 (see table 3) indicates the compliance method presented herein is 
probably as accurate as any available. 


Teeth 



TIP S/P ROOT 

Distence Along Line of Action/Pitch 

Figure 8. - Normalized beam displacement at load. 


if = 20°; A = 1.0/P; D = 1.25/P; r + 0.38/p; 
plane strain (refs. 9 and 10). 



Tip 


S/P 


Root 


Oistence Along Line of Action/Pitch 

Figure 9. - Deflection effective fillet angle versus 
load position and gear teeth. 
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Table 3 Comparison of Tooth Beam Compliances: Full Depth,/* = 4 Teeth; ^ = 20°, A = 1.0/P, D = 1.25/P,r = 

.38//*; references [9 and 18] 


Case 

Deflection Under 1000 Lbs — lO"® In. 

Normalized Compliance - 

VT EWL 

Tooth 

Load 

Chakraborty and Hunoshikatti (9) 

Chakraborty 

Premilhot 

Equations 

4&7 

Combination 

Gear/Pinion 

Position 

S/P 

(Vti + VT2' 

^ Local 

<VBi +VI:, +yB2'*''F2> 

•v'bi +v'fi) 

Hunoshikatti 

& Baronet 

30/30 

0 

431 

139 

292 

146 

4.37 

3.40 

5.25 

50/50 

0 

380 

102 

278 

139 

4.16 

2.90 

4.85 

70/70 

0 

405 

96 

309 

154 

4.61 

2.70 

4.50 

30/30 

2.5 

741 

135 

606 

- 

(Pair) 

18.2 

(Pair) 

~I3.0 

(Pair) 

18.0 


Gear Tooth Stress Sensitivity 


Derivation of Stress Formula. 


For incorporation in the High Contact Ratio Dynamic Gear Program (ref. 1) a gear tooth stress 
sensitivity or stress per load expression was needed that would be applicable for both LCRG and 
HCRG. After a review of various methods, it was decided that the Hey wood formula (ref. 4) was the 
most viable; however, it was modified so that it could be more easily applied and yet give 
approximately the same stress sensitivities as the Kelley/Pederson and Hey wood methods. 
References 13 and 4 had shown that these latter two methods correlate well with experimental and 
finite element analyses of standard gear teeth. 

The modified Heywood formula for tooth stress sensitivity is 


atV 


L Cosit>, 




CD @ 


+ 




® 0 


Tan0z, 1 

hs J 
© 


( 18 ) 


where '/4 per Heywood (ref. 4) and the rest of the nomenclature is defined in figure 10. The values 
of hs, l^, and are obtained from the gear tooth geometry, the load position, and the point of peak 
stress in the fillet defined by angle js- This formula is based on simple stress analysis and uses 
parameters that are easier to define than those used in the Heywood (ref. 4) and Kelley/Pederson 
(ref. l^formulas. Equation (18) consists of five factors, in order, stress concentration of the 
fillet, @ beam cantilever bending stress, bending load proximity stress, axial load proximity 
stress, and (?) axial stress. The formula logically makes the proximity effects die out as the bending 
and axial load get further from the fillet and closer to the neutral axis. 

In contrast to the Heywood and Kelley/Pederson methods, for the modified Heywood method 
the position of the maximum stress in the fillet, which is defined by 75 and determines the values of 
hs, Ig, and Z^, is determined by maximizing the stress neglecting the proximity terms in equation 
(18) — terms and (^. This is done by taking the derivative of equation (18) without the proximity 

terms, setting it equal to zero, and solving for 7 ^. The resulting transcendental equation is 


Tan 7 j,+ | = ^1 + A6Ai ''^Ai/ [fij ^4+ .416/li ’'^ 

- (l/3 + .016/li ’).4, Tan4>i] (19) 

where A j= ho/r + 2(1 - Cos and Bj = Ig/r +Sin 7 ^,. (See fig. 10.) Equation (19) can be easily 
solved by iteration in the given form by assuming a value for 7 ^, in the expressions for A, and Bj and 
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Figure 10. 
formula. 


Nomenclature for modified Heywood Figure 11. - Fillet angle for maximum stress versus 

geometry parameters t/hg and r/hg. 


solving for 7 j, + 1 - Usually only two or three iterations are needed to obtain 7 ^ to within 0.1 Some 
additional refinement of equation (19) should be done because it gives 7 j «0 for very large values of 
IJr, whereas it should be about 20°. For short beams, such as gear teeth, this deficiency is probably 
minor. 

Formula (19) was solved for the simple cantilever beam case with no axial load, that is, </>£ = 0. 
The position of the maximum stress as defined by angle 7 ^ is given in figure 1 1 as a function of fillet 
radius, r/hg, and load position, T/hg. The results show that for typical LCRG values of r/hg and l/hg, 
the angle of about 30° or so assumed by Heywood is a fairly reasonable average value. However, for 
HCRG with greater l/hg values, the magnitude of 75 for maximum tooth stress should be 
considerably lower than that used by Heywood. Figure 1 1 shows that the smaller the fillet the higher 
up the fillet the peak stress occurs, which is consistent with Jacobson’s test results given in reference 
14 and Winter’s test work in reference 18. In general, the point of peak stress varies significantly with 
load position and geometry and, therefore, should be Included in the derivation of the stress 
sensitivity of gear teeth. 

Evaluation of Stress Formula. 

Various gear tooth designs were evaluated using the modified Heywood formulas (18) and (19) 
as well as those developed by Heywood (ref. 4) and Kelley /Pederson (ref. 13). Comparison of the 
results with photoelastic measurements is given in figure 12. This figure shows that all three methods 
do a good job in predicting the maximum tensile tooth stressing as a function of load position, 
geometry, and number of teeth. It appears that the modified Heywood method correlates slightly 
better than the other two in some instances. Because of this and its simplicity in defining the 
necessary parameters, its use is justified over the other methods. In general, the correlation with the 
photoelastic results are within about 5 percent, and slightly on the conservative side. 
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Figure 12. - Improved modified Heywood results versus photoelastic tests, 
Kelley/Pederson and Heywood. 


Table 4 Comparison of Improved Modified Heywood and Exact Methods for Full Depth Involute Spur Gear Teeth: /4 = 
1/P, D = 1.2S/P, r^. = .380//», r^- = .433//»; reference [15] + 


Ca» 

Pitch Line Loading 
- 0.0/P 

Tip Loading 
H„ = 1.0/P 



Stress Factor 

Peak Stress 

Stress Factor 

Peak Stress 
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Jt 

+ 
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wM 

of 

"Exact" 
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• 
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• 

Heywood 

• 

Heywood 

• 

Heywood 


20 

0.328 


1.17 

_ 

0.221 


1.14 

- 


30 

- 


- 

1.11 

- 


- 

1.01 

20° 

40 

0.407 



1.11 



1.11 

1.01 


60 

0.462 


bh 




1.11 

1.01 


80 

0.473 


1.15 




1.11 

1.01 


150 

0.483 


1.14 

1.13 



1.09 

1.00 


20 

0.451 

_ 

1.24 

_ 

0.307 

0.254 

1.19 



30 

- 

0.498 

- 

1.13 

- 


- 


2S° 

40 

0.494 

0.531 

1.19 

1.13 

0.315 

0.290 

1.15 



60 

0.517 

0.565 

1.17 

1.13 

0.317 

0.298 

1.14 



80 


0.586 

1.17 

1.14 


0.304 

1.14 



150 

0.569 

0.608 

1.17 

1.14 


0.307 

1.14 

0.99 1 


CJ| • L| P/J|W; W ■ Width? P » Pitch; ■ Tangent Load at Pitch Radius 


Table 4 and figure 13 present the comparison of the stress sensitivities for a family of 20° and 
25° pressure angle, full depth involute gears using the modified Heywood and an “exact” 
transformation analysis (ref. 15). In general, the trends with number of teeth, load position, and 
pressure angle are predicted by the modified Heywood formulas. For the 20° pressure angle gears, 
the pitch circle loading results correlate very well, but the tip loading results are about 10 percent 
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Table 5 Finite Element, Modified Lewis, and Modified Heywood Tooth Stressing: Rp = 2.50', 0/. = 20°, W 

1.0", A = .25', D = .289', r = .060', = .390' 


Par«m«t«r 
S— Figur* 10 

Finita Eiamant 
& tmpuiss Tast 

Modifiad Lawis 

Modifiad Haywood 

Load Position 

1 

2 

3 

1 

2 

3 

1 

2 

3 

— in. 

~ 2.70 

~2.55 

^2.40 

~2.70 

'2.55 

'2.40 

2.70 

2.55 

2.40 

0’l 

27.10 

18.83 

6.62 

27.10 

18.83 

6.62 

27.15 

18.88 

6.68 

hj - in. 

- 

- 

- 

0.45 

0.46 

0.47 

0.426 

0.434 

0.451 

J?s 

- 

- 

- 

0.40 

0.26 

0.15 

0.403 

0.255 

0.150 

0,/L - Pfi/lb 

16* 

12.8* 

10.6* 

15.3 

12.2 

9.3 

20.2 

14.7 

11.7 

(7s + r»-d.* 

^10* 

~30* 

^45* 

- 

- 

- 

23.0 

32.5 

47.1 


Faiiur* Test; '*'Be(ow Surface 0.008"; Therefore Values Low 



Figure 13. - Comparison of geometry factors from 
various formulas and "exact" analysis. <!’ p = 

20°; P = 1; A = 1/P; D = 1.25/P; r = 0.38/P; a = 
LtP/JtW. 


high, whereas the reverse is true for the 25° pressure angle gears. The “exact” analysis indicates the 
peak stress occurs further up the fillet than the modified Heywood analysis results. (See table 4.) The 
accuracy of the analytical function in defining the fillet region could influence its results and the 
correlation with the modified Heywood results. 

Wallace and Seireg performed a finite element analysis of a gear tooth (ref. 16) and also 
conducted impulse strain gage tests of a standard, 4 pitch gear tooth. Comparison of their finite 
element results, those by the modified Lewis formula, and those by formulas (18) and (19) are 
summarized in table 5. This summary indicates that the modified Lewis results are a little low, 
particularly if the finite element results are corrected for the effect of element size. In contrast, the 
modified Heywood results are slightly higher, 10 to 20 percent, but this difference would be 
decreased if the effects of the coarse finite element mesh were taken into account. The locations of 
the peak stressing predicted by the modified Heywood method for the three tested tooth loading 
positions correlate well with those predicted by the modified Lewis method and the failure locations 
found during the impulse tests. Some of the discrepancy in the results between the two methods of 
analysis appears to be due to the definition of the tooth geometry. Although the distance from the 
load to peak stress, 1^, is about the same for the two methods, their respective tooth thicknesses, h^, 
are about 5 percent different, which would account for about 10 percent of the differences in the 
stress sensitivies. 
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As a more fundamental check on the modified Heywood stress sensitivity analysis, it was used to 
calculate the maximum stress and its position for photoelastic models H and F tested by Dolan and 
Broghamer (ref. 17). The first model is a simple rectangular, cantilever beam loaded normally close 
to the fillet, J/ho = 0.55, whereas the second model is similar to the first but the beam is tapered 15° 
on each side with the load normal so as to introduce an axial compressive load. The photoelastic 
results for the first model gave a tensile stress sensitivity of aW/L = 6.2Q/'\n. and a maximum stress 
position at about 25° up the fillet. The modified Heywood method gives a stress sensitivity value of 
6.06/in., 2 percent low, at an angle of 24° up the fillet. The photoelastic results for the second model 
gave a tensile stress sensitivity of 4.46/in. at about 35° up the fillet, whereas the modifed Heywood 
method predicts values of 4.37/in. and 33°, respectively. The good correlation in magnitude and 
position of the maximum tensile stress shows that formulas (18) and (19) can handle accurately loads 
close to the fillet. 

The importance of having an accurate method for calculating the maximum stress in various 
gear teeth was illustrated in reference 1. A comparison was made of the calculated tooth stress 
sensitivity versus load position for both HCRG and LCRG using various conventional and Heywood 
type analyses. The results, reproduced in figure 14, show that the three Heywood type analyses agree 
quite well with each other, whereas the usual Lewis and AGMA analyses give lower and significantly 
different results. The relative differences between the modified Heywood and the Lewis and AGMA 
results at the load position for maximum likely peak stress are quite different for LCRG and HCRG, 
being about 20 to 28 percent higher for the LCRG and only about 10 to 18 percent higher for the 
HCRG, respectively. Such differences would have significant effect on the tooth bending fatigue life. 


Improvements in Dynamic Gear Program 

The High Contact Ratio Dynamic Gear Program developed and discussed in reference 1 already 
incorporated the stress sensitivity formula developed herein. However, the program since has been 
improved by including a preprocessor which defines the gear tooth structural and manufacturing 
parameters from a few simple input parameters and calculates the gear tooth pair compliance based 
on the equations developed herein. A postprocessor has been added which calculates the design 
outputs of PV, Hertz stress, flash temperature, load, and stressing for the two gears. In addition, the 
formulation of the basic program was expanded to include centerline error, which affects the 
pressure angle and length of engagement, and a three term definition of tooth profile relief. A, rather 
than a two term expression (see ref. 1), for example, for engagement 

A, = C,(S^-S,)2[l+y,(S^-S,)l (20) 

in place of 
a, = C,(S„,-S,)2 


Effect @Peak Stress Position 
•^mod Hey^'^Lewis “1-28 
“mod Hey/“AGMA“’-2° 

^eyllVOod 



®oe ^ SqjI 

Position On Line Of Action 
a: 8 Pitch/32 Tooth LCRG 



^oe ® 

Position On Line Of Action 
b; 8 Pitch/32 Tooth HCRG 


Equation 20 



Figure 14. - Comparison of gear tooth stress 
sensitivities for different analysis methods. 


Figure 15. - Tooth profile modification shapes; 
engagement cam. 
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Figure 17. - Dynamic tooth loads and stressing for HCRG-tooth error case: 2.2 
contact ratio, 8 pitch, 20° pressure angle; 5000 in-lb torque; 44/52 gear 
ratio; 6-in. center distance; ref. 1 analysis. 


The effect of this additional parameter yg on the allowable profile or cam shapes is shown in figure 
15. Figures 16 and 17 are examples of the plotted outputs for the improved program. The first figure 
depicts the results for an HCRG, 8 pitch, 44/52 tooth pair with a nominal 2.20 contact ratio, a 20° 
pressure angle, and no tooth error under 5000 in. -lb. of torque. The second figure is for the same gear 
pair but with tooth spacing error and is for eight mesh cycles. 

Conclusions 

Formulas have been developed based on work done by Heywood (ref. 4), Weber (ref. 2) and 
O’Donnell (ref. 3) for the stress sensitivity and compliance of low and high contact ratio, involute, 
spur gear teeth. The parameters used in these formulas require the derivation of the effective fillet 
length or angle rather than assuming particular values as done by Heywood (ref. 4) and O’Donnell 
(ref. 3). The stress sensitivity formula is a modified version of the Heywood formula, using simple 
beam parameters. The compliance formula uses O’Donnell’s foundation flexibility factors and 
Weber’s local contact compliance. 

The modified Heywood gear tooth stress sensitivity formula was evaluated by comparing its 
predicted results with test, finite element, and analytic transformation results for gears with various 
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pressure angles, tooth proportions, number of teeth, and load contact position. In general, the 
formula predicts the maximum tensile fillet stress within about 5 percent and gives the proper trends 
with contact position, number of teeth, and tooth geometry. It also predicts fairly well the location of 
the peak stress in the fillet. Comparison with similar formulas by Kelley/Pederson and Heywood 
show good correlation, but with the proposed modified Heywood method giving slightly more 
consistent results as well as being easier to use. 

The beam and foundation portion of the gear tooth compliance formula was evaluated using 
available test, finite element, and analytic transformation results and found to correlate quite well. 
The inclusion of a derived, variable effective fillet length or angle appears to improve upon the results 
which assumed a given effective fillet length. In general, the combined beam and foundation 
deflection as a function of gear teeth and load contact position agreed well with other methods. The 
Weber method of calculating the local contact deflection is universally used and was found to give 
more consistent results than two other methods. The gear tooth pair compliance formula was found 
to predict the overall compliance (beam + foundation + local) consistent with the test results of 
Howland. 

In summary, although further evaluation of the gear tooth stress sensitivity and compliance 
formulas presented herein might be desirable, the evaluation results indicate that the formulas predict 
the spur gear tooth pair stress sensitivity and compliances quite well. 
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Appendix — Definition of Different Gear Tooth Forms 


General equations are given which define the involute tooth form based on pitch, P\ pitch 
pressure angle 4>p\ pitch radius, Rp\ circular pitch tooth thickness tp\ root radius, R/, fillet radius, r; 
and tooth width, W. The load position on the tooth can be defined in terms of the distance along the 
line of action from the pitch circle, S, or by the contact or load radius Rp (fig. 2). Knowing these eight 
quantities, the basic tooth parameters common to all three tooth forms (fig. 1 ), are 

number of teeth: N=2PRp 

base radius: Rg = Rp Cos 4>p 

base radius included tooth angle: aB = tp/2Rp + lan 4>p-<i>p 

contact radius: R^- S^ + 2Rp Sin 4>pS + Rp 

contact pressure angle: Cos <j>L = RB/RL 

contact included tooth angle: a/^ = « 5 -Tan + 

contact tooth thickness: hi = 2Ri Sin 

contact load angularity with L: <j)[ = {(fti - a^) 

Gear Tooth Form S; standard case: R^^+ 2rR^>R\ 

W2 = (R, + r)2 - 2rV (Rr + r)2 -R^ + r^ 

Cos ^ = Rb/R 

a = aB — Tan 4> + 4> 

h = 2R Sin a 

y = <t)'={^-a) 

Io = [Rl Cos ai-{hi Tan (p[)/2-{R Cos a + r Sin 7 )] (stress) 
ho = 2R Sin a - 2r(l - Cos 7 ) = h - 2r(l - Cos 7 ) 

The term 7 ^ is iterated using equation (19). 
hs = ho + 2r(l - Cos 7 ^) 
l's=lo + r Sin 7 ^ 

/, = /; + (/! Tan 4>'lV2 

The term 7/7 is iterated using equation (7) as in. 

X = R Cos a-R Cos a 

J=Ri Cos ai~{hi Tan 4>[)/2-R Cos a (compliance) 
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The term /?£ is iterated until Cos a[-R Cos a) = /. 
hl = 2Rl Sin a£ 

Gear Tooth Form A; tangent Case: R^+2rRr<Rg (fig. 18(a)) 

Sin ar=r/Rf + r) 

hg = 2(i?;. + /-)Sin + a^) - 2r 

Iq = Rl Cos a^-ihi Tan <f>'j)/2-{Rr + r) Cos (a^ + a;.) 

/zg = 2R B Sin ( ~ "y) = 2/?g Sin ot g 

Between the tangency point of the fillet and base circle, for example, R<R<Rb, 

h = [(^r + r)/Cos a, + X/Cos 7| 

Rg L J 

g 

= hg , where A" = (R — R) Cosag 

Rb 

ag is iterated using equation (19), and hg, Ig, and, Ig are the same as for form S 
agis iterated using equation (7), but including additional beam section (/?g-/?) Cos ag 
Gear Tooth form B, parallel case: R^^+2rRr<R ^ b (f^8. 18(b)) 

7= —6 

e=-K/Nov 180V7V 
8 = Rb Sin (0-ag)-/? 

A = Rb Cos (0-ag) -(/?;. + /■) 

Iq = Rl Cos u^-Rb Cos ag + A Cos 6 + r Sin 6 — (hi Tan <A£)/2 
ho = hB-1- A Sin 6-2r(l -Cos 6) 
hB = 2,RB Sin ag 

Between the tangency point of fillet and base circle, for example, 0<A'< A Cos 6, 

jg is iterated using equation (19), and hg, Ig, and Ig are the same as for form S. 

7 gis Iterated using equation (17) but including beam section A Cos 6 
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